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Abstract 

A conservative Newton system q = —'W{q) in is called separa- 
ble when the Hamilton-Jacobi equation for the natural Hamiltonian 
H = ^p^ + V{q) can be solved through separation of variables in some 
curvilinear coordinates. If these coordinates are orhogonal, the New- 
ton system admits n first integrals, which all have separable Stackel 
form with quadratic dependence on p. 

We study here separability of the more general class of Newton 
systems q = — (cof G)~^'VW{q) that admit n quadratic first integrals. 
We prove that a related system with the same integrals can be trans- 
formed through a non-canonical transformation into a Stackel separa- 
ble Hamiltonian system and solved by qudratures, providing a solution 
to the original system. 

The separation coordinates, which are defined as characteristic 
roots of a linear pencil G — of elliptic coordinates matrices, gener- 
alize the well known elliptic and parabolic coordinates. Examples of 
such new coordinates in two and three dimensions are given. 

These results extend, in a new direction, the classical separability 
theory for natural Hamiltonians developed in the works of Jacobi, 
Liouville, Stackel, Levi-Civita, Eisenhart, Benenti, Kalnins and Miller. 
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1 Introduction 



The method of separation of variables for natural Hamiltonian systems is 
well understood [HE]- In terms of (orthogonal) separation coordinates, the 
natural Hamiltonian H{q,p) = + V{q) becomes 

n 

(1.1) H{x,y) = H{q{x),p{x,y)) = \Y, H'^x) + V {q{x)) , 

i=l 

and separation means that its Hamilton-Jacobi equation, obtained by substi- 
tuting yi = dS{x; a)/dxi, admits an additively separated solution S{x; a) = 
J2^=i ^ii^i'^ After this substitution the Hamilton-Jacobi equation splits 
into a system of first order ODEs for the functions Si^xf, a) and the solution 
can be expressed by qudratures. 

The well known Stackel theorem ^Hl ^] gives a necessary and sufficient 
conditions for H^^, . . . , and V{q{x)) to admit a separated solution of the 
Hamilton-Jacobi equation. This Stackel condition is purely algebraic, and we 
say that the Hamiltonian ()1.1|) has separable Stackel form if these conditions 
are fulfilled. The change of coordinates q{x) to separation coordinates x is 
a priori not known and the problem of finding these variables for any given 
potential V{q) has been solved only recently |14j . 

For general dynamical systems x = f{x), x = {xi, . . . ,x„)"'", or for non- 
canonical Hamiltonian systems z = n{z)VH{z), where n{z) is a Poisson 
matrix, the notion of Hamilton-Jacobi separability is not well understood. 
The reason is that for non-canonical Hamiltonian systems there is no natu- 
rally associated Hamilton-Jacobi equation that can be solved by an additive 
ansatz S = ^"^^ Si. 

In this paper we study systems of autonomous second order ordinary 
differential equations of the form q = M{q), q = (gi, . . . , qn)'^ , where accel- 
eration q equals a force M{q), which does not depend on velocities q. Such 
systems have been given the name Newton systems. We consider Newton 
systems of a special type having a force M(g) that is quasi-potential in two 
different ways M{q) = - {cof G)-^VW{q) = -{coi G)-^VW{q), where G, 
G are elliptic coordinates matrices, cof G = (detG)G~^ and W{q), W{q) 
are quasi-potentials. Such cofactor pair systems naturally generalize separa- 
ble potential Newton systems q = —W{q) in i?" since the function W{q) 
becomes an ordinary potential for G = I. 

Any cofactor pair system can be written as a dynamical system q = p, p = 
M{q), and is known to possess, in the extended phase space, a non-canonical 
Hamiltonian formulation z = n{z)'VH{z), z = (gi, . . . , g„,pi, . . . . 
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These Newton systems have been proved to be integrable in a somewhat non- 
standard way through embedding into a Liouville integrable system |7l[TT|. 

In the first part of this paper we show that any cofactor pair system 
q = — (cof G)~^'VW{q) = — (cof G)~^VW (q) has a related system that can be 
transformed through a non-canonical transformation into a canonical Hamil- 
tonian system 

dH(x,y) dHix,y) 

^ = S ' y = 5 ' 

oy ox 

which has separable Stackel form fTheorem 12.11) . A solution of the related 
system then provides a solution of the cofactor pair system by rescaling the 
time variable. This result indicates a simple and natural way of prescribing 
the Hamilton- Jacobi separability property to dynamical systems that are not 
canonical Hamiltonian systems. Such systems can be considered to be separa- 
ble whenever there exists a transformation (usually not a canonical one) into 
a canonical Hamiltonian system that admits the classical Hamilton-Jacobi 
separability (see Remark 12.81 and Remark 12. This observation allows for 
extending the concept of separability and the concept of solving through 
separation of variables to large classes of dynamical systems that otherwise 
could not be considered separable in the Hamilton-Jacobi sense. 

Separability for generic cofactor pair systems has already been studied 
and separation coordinates in the extended phase space of variables 
(g, p, d) have been derived using the theory of bi-Hamiltonian systems devel- 
oped by the Milano school 0. Separability for the more degenerate class 
of driven cofactor systems has also been studied [H]. In this paper we ap- 
proach the problem of separability without using the concept of extended 
phase space and find a direct transformation into separation coordinates in 
which the related system attains separable Stackel form. The advantage of 
this approach is that it explains, on the level of differential equations, the 
mechanism of separability and gives simple tools for solving these systems 
without resorting to the geometrical language of the Milano school. 

The second part of this paper studies cofactor-elliptic coordinates defined 
as characteristic roots of a linear pencil G — /iG of two elliptic coordinates 
matrices G, G. They separate cofactor pair systems and appear to be a 
natural, however more complicated, generalisation of the elliptic coordinates. 

Elliptic coordinates Xi, . . . , x„ play a central role in the classical sepa- 
rability theory for natural Hamiltonians (jl.lj) . They are defined as zeroes 
of 

i+E7zv=n(--^.)/n(^-^^) 

i=i * j=i I k=i 

where Ai, . . . , A„ are parameters. The elliptic coordinates are mother of 
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all orthogonal separable coordinates on for the natural Hamilton- Jacobi 
equations. This means that all other separation coordinates can be recovered 
from the elliptic coordinates through certain proper or improper degenera- 
tions of the parameters Aj. For cofactor-elliptic coordinates we find a formula 
[see ()4.3j) ]. which closely reminds of the above formula but encompasses con- 
siderably more coordinate systems. These coordinates are usually not orthog- 
onal in the ordinary Euclidean sense, but instead orthogonal with respect to 
the scalar product (f , w) ^ v^Gw. For instance, in two dimensions, cofactor 
eUiptic coordinates xi, X2 are defined as solutions of 

(gi + ^xf ^ ql ^ ^ 

Ai — (1 — e)x + e'^x^ A2 — x 

For e = this equation defines two-dimensional elliptic coordinates. Thus 
Ai, A2 play a similar role as before, while for other values of the parameter e, 
the curves of constant value of x cover the plane in a complicated way. Their 
pattern strongly depends on the relative values of Ai, A2, and usually not 
all points of the plane are parametrised through cofactor elliptic coordinates 
Xi, X2- In such domains separation takes place for complex values of Xi, X2- 
Several examples of these coordinates are given and illustrated by pictures. 

In this paper we discuss and exemplify only the generic type cofactor- 
elliptic and cofactor-parabolic coordinates that play a fundamental role for 
cofactor pair systems. This study is the first step into a new fascinating world 
of nonorthogonal separation coordinates defined by families of non-confocal 
quadrics. 

1.1 Properties of cofactor pair systems 

Before discussing separability of cofactor pair systems, we need to define 
them and recall some facts from [7j. 

An elliptic coordinates matrix is a matrix-valued function of g G i?" that 
can be written in the form 

G{q) = aqq^ + qjS^ + I3q^ + 1 

[we refer to elements q = (gi, . . . , qnY of -R" as vectors and view them as n x 1 
matrices], where a G -R is a constant scalar, j3 G i?" is a constant vector and 
7 = 7""" is a real constant symmetric nxn matrix. By construction, G is also 
symmetric. 

A simple but very useful fact is that an arbitrary linear combination 
of elliptic coordinates matrices is an elliptic coordinates matrix too. This 
implies that all statements to be made about elliptic coordinates matrices 
also hold for pencils G — fiG of elliptic coordinates matrices. 
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For a quadratic matrix A, we write cof A to denote the corresponding 
cofactor matrix, which is defined by (cof A)^ = A{coi A) = {detA)I. 
To a given eUiptic coordinates matrix G, there is an associated vector 

N = aq + (3= IVtrG 

[V denotes the gradient operator dg = {d/dqi, . . . ,d/dqn)'^ expressed in 
Cartesian coordinates]. If G is non-singular, it is possible to represent this 
vector as 

(1.2) Ar = iGVlogdetG 

[we use logarithmic derivatives like VlogF formally to denote (1/F)VF 
regardless of the sign of F] . The associated vector satisfies 

(1.3) V (A^"^ (cof G')iV) = 2a (cof G')A^. 

A Newton system q = M{q) in i?" is called a cofactor pair system if the 
force can be generated in two different ways as 

M(q) = -(coiG)-^VW = -(coiCy^VW, 

where 

G{q) = aqq" + q(3' + (3q' + 7, det G ^ 0, 

G{q) = aqq^ + qp^ + pq^ + 7, det G ^ 0, 

are two linearly independent nonsingular elliptic coordinates matrices. Here 
W and W are two functions on i?" called quasi-potentials. Clearly, these 
quasi-potentials have to satisfy 

(1.4) (cof G)-^VW = (cof G)-^VW. 

The Frobenius compatibility conditions for this system of equations, rewrit- 
ten in terms ol K — W/ detG and K — W / det (7, are referred to as the 
fundamental equations. 
By defining 

E =\q'^ {cof G)q + W and £^ = |g"^(cof G)^ + VF, 

two quadratic first integrals of "energy type" are constructed. In fact, the 
so-called "2 ^ n theorem" states that any cofactor pair system admits n 
quadratic first integrals 

E^^^ = ig^A^^'^g + ly(^), k = 0,...,n-l, 
6 



where the matrices A^''^ are given by the generating function 

n-l 

coi{G + i^G) = Y,!^''^^'"^- 

k=0 

The quasi-potentials W^''^ are given as solutions to the differential equations 

which are integrable if the fundamental equations are satisfied by K — 
W/ det G and X = W/ det G. In particular, E^^^ = E and = E. 

In this paper, we focus on cofactor pair systems having the property that 
the polynomial det(G' — //G) has n functionally independent roots fi = Uk{q). 
We call such systems generic. In the generic situation, the functions Uk{q) 
are necessarily non-constant and give rise to the non-singular Jacobian 

J= (Viii,..., Vit„), det J 7^0. 

It is thus clear that we can define new coordinates as being the roots = 
Uk{q). For this change of variables we can write the chain rule compactly as 
V = Jdu- The gradients Vuk are G-orthogonal, (V%)^GVMfc = for j ^ k, 
and give rise to generalized metric coefficients Ak — {'Vuk)'^G'Vuk- Thus we 
have the generalized orthogonality relation 

(1.5) J"^G'J = ^ = diag(^i,...,Z\„). 
Writing 

X = G-^G, 

we see from G — fiG = G{X — /.jJ) that are eigenvalues of X. The 
corresponding eigenvectors arc the gradients Vwfc, so that XVuk = u^Vuk 
for all which we can write as 

(1.6) XJ = JU, U = diag(Mi, . . . , m„). 
This implies that X satisfies 

(1.7) XVlogdetX = VtrX 
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2 Main separability theorem 



In jZI it was shown that to any cofactor pair system q = M{q) there is a 
related bi-Hamiltonian system 

iG p \ fdA 

-iG"^ l{Np^-pN^) M + dN \ \dp\ (ddetG) 
V -{M + dN^ / W 

in the extended {2n + l)-dimensional phase space obtained by taking p = q 
and introducing an extra variable d. On the hyperplane d = 0, the bi- 
Hamiltonian system reduces to 







p 


H 


\d) 







(2.1) 

This system has the same trajectories as the system 
(2.2) 





which clearly is equivalent to q = M{q) when d = 0. Since the trajectories 
are the same, the solutions to ()2.1|) and ()2.2|) differ only by a scaling of the 
independent variable. So if the solution of ()2.1|) is known, the solution of 
()2.2j) can be obtained by substituting the correct expression r(t) for r. In 
this sense, every cofactor pair system is equivalent to a Hamiltonian system. 

Using these ideas, it has been shown that any cofactor pair system is 
Liouville integrable (in a somewhat non-standard sense), provided that the 
first integrals in the 2 ^ n theorem are functionally independent and in 
involution. 

For cofactor pair systems, explicit integration through the Hamilton- 
Jacobi method has also been shown TD^ . Guided by the theory of Darboux- 
Nijenhuis coordinates for bi-Hamiltonian systems, the separation coordinates 
{u, s, c) were introduced, where Uk{q) are defined as roots of the polynomial 
det(G + /iG), while 

and c = (det G)d. In these coordinates, the related system ()2.1|1 takes the 
form 
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and the Hamiltonian h{u, s, c) attains separable Stackel form. When c = 
one then obtains a 2?T,-dimensional canonical Hamiltonian system in {u, s) 
variables, so the classical theory is applicable. The connection between the 
(2n+l)-dimensional system and the cofactor pair system q = M{q) is however 
intricate and it is not easy to see the separation mechanism for cofactor pair 
systems. 

We will now give an alternative formulation of this result, in the spirit of 
the classical Stackel approach. The techniques used in the proof are further 
extensions of the results in [Hj. 

Theorem 2.1 (Stackel separability for generic cofactor pair sys- 
tems). Let (Pq/de = M{q) = -{cof G)-^VW = -{coi G)-^VW be a cofac- 
tor pair system in i?" written as a dynamical system 



(2.4) 

with the related system 



d 



— fA = (p 

dt \pj \M 



obtained by rescaling the vector field ()2.4j) . 

Suppose that the polynomial det(G' — fiG) has n functionally independent 
roots fi = Uk{q), so that they define new configuration coordinates Uk = Uk{q)- 
Define new momenta coordinates as 

(2.6) ^'^ = ^'=(^'P) = i;E^P^ 

where Ak = {VukY GVuk are generalized metric coefficients. In these new 
coordinates we have: 



1. The related system ()2.5|) has a canonical Hamiltonian formulation 
with Hamiltonian 

n 

(2.8) H =\{deiG)^Aisl + W{q{u)). 

1=1 
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2. The Hamiltonian ()2.8|) has separable Stdckel form 

1=1 ^ 

where fi and gi are some functions of one variable. The Stdckel matrix 
is the Vandermonde matrix ()A.6|1 and consequently, the coefficients 
Co{ui)/U'{ui) are elements of its inverse (see Appendix IX)) . 

3. Systems ()2.4|) and ()2.5|) have a common set ofn quadratic first integrals 

^(.) ^ ^ ^n_^-iM ,2 ^ k = 0,...,n-l, 

^ U'[ui) 

where fi and gi are the same functions as above. These first integrals 
are functionally independent and in involution. The Hamiltonian ()2.8|1 
is given by H = E'^'^~^\ 

4. The Hamilton- J acobi equations E^''\u,duS) = at, k = 0, . . . , n — 1, 
admit a common separated solution S{u;a) = Sk{uk', o^) . In partic- 
ular, S is a solution to the last equation 

i(detG')5^Z\J— j +W{q{u))=an-i 
1=1 ^ * ^ 

associated to the Hamiltonian fj2.8|) . The corresponding Hamiltonian 
system ()2.7j) and the cof actor pair system ()2.4|1 can therefore be solved 
by quadratures. 

We shall clarify the content of this theorem and the connection with clas- 
sical separability of natural Hamiltonian systems in the subsequent remarks. 

Remark 2.2 (Change of coordinates). It is convenient to rewrite for- 
mula ()2.6|) for the change of momenta coordinates in matrix notation. With 
the previously introduced generalized metric matrix A and the Jacobian J, 
we have 

s = A-^J^p. 

It is easy to see that if J is invertible, then so is the full phase space trans- 
formation {q,p) {u,s). 

Remark 2.3 (Definition of momenta). One can show that definition ()2.6|) 
agrees with ()2.3|) . The form ()2.6|) is however more suitable for our proof. 
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Remark 2.4 (Comparison with the classical case). Any cofactor system 
q = — (cof G)^^ Viy admits a first integral of energy type, E = |g"'"(cof G)q + 
W, by construction. Stated in a non-symmetric way, Tlieorem 12.11 savs tliat 
a cofactor system q = —(cof G)^^'VW is separable if there is an extra first 
integral, also of energy type, E = ig"'"(cof G)q + W such that the eigenvalues 
of G~^G are functionally independent. 

A potential Newton system q = —W is of cofactor type with G = I, 
and the accompanying first integral is the total energy E = \q^ q + V . The 
theorem therefore says that it is separable if there is an extra first integral E 
of cofactor type with functionally independent eigenvalues. Its existence also 
implies the existence of a set of n first integrals. 

Benenti [T] has shown that the separability of a potential Newton system 
is equivalent to the existence of a very special first integral, corresponding 
to a "characteristic" Killing tensor K. For instance, in the case of systems 
separable in elliptic coordinates, this tensor is expressed through the "inertia" 
tensor G = —qq'^ + diag(Ai, . . . , A„) as = (trG)/ — G. It has also been 
noticed that the eigenvalues of G are separation coordinates in this case. 

In our terminology, G is an elliptic coordinates matrix, and one can show 
that the eigenvalues of G define elliptic coordinates (see Proposition 14. 2p . 
which also justifies the terminology. Moreover, by the construction in the 
2 ^ n theorem Theorem 5.3], A^^-^) = (trG)/ - G. When n = 2, we 
thus have K = A^^^ so that Benenti's characteristic first integral is the same 
as the first integral E considered here. However, when n > 2, the two types 
of first integrals are different, even though they both imply separability for 
the classical potential Newton system. 

Remark 2.5 (Novel features). By using generalized orthogonality ()1.5j) . 
the transformation formula ()2.(ij) for the momenta coordinates can be written 
as 

s = (Gjy^p. 

This can be compared with the corresponding formula for a canonical trans- 
formation of coordinates: a change of configuration coordinates g m in- 
duces a change of momenta coordinates p —>■ s according to 

s = J~^p. 

Thus we see that the change of coordinates considered in this paper is non- 
canonical, unless G = I, which, as mentioned before, specializes the theory 
to the classical case of separable potential systems. 

Another novel feature of this separability theorem is the fact that a related 
system, having the same trajectories, is solved instead of the original system. 
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Again, if G = /, the two systems coincide, and the theory speciahzes to the 
classical case. 

Remark 2.6 (t h-^ r). Solutions of the two systems ()2.4j) and ()2.5|) are related 
by the function r(t) defined as a solution to the differential equation 

dr 

(2.9) detG(gor)- — = 1, 

where g is a solution to ()2.5j) . Indeed, if {q,p) solves ()2.5j) . then (g o r, p o r) 
solves ()2.4|1 as is easily verified by the chain rule: 

|(gor) = (gor)| = detG(gor).(por).^=por, 

and similarly for d{poT)/dt = M(gor). The solution of ()2.9p can be obtained 
by one further quadrature. 

Remark 2.7 {2 ^ n theorem). The quadratic first integrals referred to 
in Theorem 12.11 coincide with those given by the 2 ^ n theorem Theo- 
rem 5.3], as follows from the relation 



, -p , n—l n 

(2.10) (J) cof(G + /.G)(J) =(detG)^/.^^a._._i(^,)A-.| 

k=0 j=l 



proved below. This relation shows that the first integrals considered in the 
two theorems have the same kinetic part and therefore have to be the same, 
since the Hamiltonian ()2.8|1 has a non-trivial quadratic dependence on all 
momenta variables. 

The separability theorem thus implies the 2 ^ n theorem in the case of 
generic cofactor pair systems (since the separation coordinates u are possible 
to define for such systems). See further Section lOl 

Proof. We have that dq/dt = p, if t and r are related as in Remark 12.61 Thus 

dt 

Further, the eigenvalue-eigenvector relation ()l.fi|) gives GJ = GJU, so that 
J^GJ = Z\ZY by generalized orthogonality ()1.5|) . It is now clear that (when 
— /i is not an eigenvalue of X) 

(^)^cof(G + ^G)(^) = det{G + fxG)s'^A{j^{G + fiG)J)~^As = 

(det G) det(X + /iJ) s"^Z\(W + /i/)"^s. 
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which can be rewritten as a sum of products 



n 



n 



(detG')5^n 



Ui + IJ, 




j=i i=i 



By using ()A.2j) . we can expand these products to find the right hand side 



Remark 2.8 (Hamiltonian separability for non-Hamiltonian sys- 
tems) . The example of cofactor pair systems indicates how to give meaning 
to the concept of separability for general Newton systems q = M{q) or more 
general second order dynamical systems. Any such system can be considered 
to be separable whenever there exist a transformation to new coordinates 
{u, s) in which the transformed equations acquire the canonical Hamiltonian 
form and the related Hamilton- Jacobi equation can be solved through a sepa- 
rated ansatz 5* = ^ Sk{uk)- Such separation can be even more general when 
the Hamiltonian for the transformed equations does not have the Stackel 
form as is necessarily the case for cofactor pair systems. 

Remark 2.9 (Separability of non-canonical Hamiltonian systems). 

The related system has a non-canonical Hamiltonian formulation in the 
{2n + l)-dimensional phase space with coordinates {q,p,d), but it has no 
directly related Hamilton-Jacobi equation. The transformation (g, p, d) —>■ 
{u, s, c) involving ()2.H|1 given in ^0] is non-canonical and gives the Hamil- 
tonian system a canonical form in {2n + l)-dimensional phase space, with 
a Stackel separable Hamiltonian. In this sense, a non-canonical Hamilto- 
nian system can be considered separable. The important thing is, that a 
non-canonical Hamiltonian system has been turned (through a non-canonical 
transformation) into a canonical Hamiltonian system for which the concept 
of Hamilton-Jacobi equations and the concept of separability is well defined. 

3 Proof of the separability theorem 

In this section, we prove Theorem 12. II The proof is rather long and technical, 
but have been structured into subsections, to help the reader in keeping track 
of the logic. To simplify the notation, we now use a dot for d/dr. 

3.1 Hamiltonian formulation 

We begin by establishing a canonical Hamiltonian formulation of ()2.5|) . First, 
we transform the equations using the transformation {q,p) —>■ {u, s). Second, 



of (PTTHI . 



□ 
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we show that the transformed equations can be expressed through the Hamil- 
tonian ^T^ . 

In order to transform the equation q = (det G)p, we use the chain rule in 
the form q = J~'^u and the definition of the momenta coordinates written as 
p = J~'^ As. When inserted, this gives immediately the equivalent form 

ii = (detG)Z\s 

by canceling J^^ . To give this equation a Hamiltonian formulation ii = dsH, 
we need to find a function H = H{u, s) such that dgH = (det G)As. This is 
an overdetermined system of PDEs, which can easily be integrated to give 

n 

H = l{detG)J2^ksl + F{u), 

k=l 

where F is an arbitrary function of u. 

We now turn to the equation p = (det G)M, and we use the representation 
M = -(cof G')-^VW^ so that 

(3.1) p = -GVW. 

The definition of the momenta coordinates gives dpSi = Z\~^VMj, so that 
p.ip is equivalent to the system 

{dps.yp = -Aj\Vu,yGVW, j = l,...,n. 

To get further, we need the following variant of the chain rule. 

Lemma 3.1 (Chain rule). 

d 



(VmjO^GV = A 



Proof. The chain rule V = Jdu together with GJ = J '^A from ()1.5j) gives 

(Vmj)^GV = {J-^Vu^yAdu, 

which equals Aj d/duj since J^^Vuj forms the j:th standard basis vector. □ 
From the usual chain rule we have 

sj = {Vsj^q + {dpSjYp, 
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which together with Lemma 13.11 shows that ()3.1|) is equivalent to 



dW . ^ 
Sj = (Vs,) 9- j = l,...,n. 

These equations can be given a Hamihonian formulation s = ~duH with the 
same Hamiltonian as above, if {VsjYq = —{d/duj){H — F) and F = W. 
The latter condition can be taken as a definition since F is arbitrary, but it 
remains to show the former condition. This is more complicated, and will be 
carried out in the next section. 



3.2 Kinetic part 

We shall show that 

(3.2) (V.,)^g = (l(det G) g A,sl^ . 

This will be done using q = (det G)p, which can be written as 

q = (det G)GJs, 
since it is possible to express p = J^'^As through 

p = GJs 

as follows from the generalized orthogonality 
We expand 

(3.3) (V^.r^=(detG)5:|jG.e|^.,. 

The derivatives of Sj can be calculated from s = A~^J'^p as 
— = — (AT^^)Gbd—Sc 



b,cA 



where the expression GJs has been substituted for p after differentiation. To 
simplify this, we use Leibniz' rule and the generalized orthogonality ()1.5|) to 
find 
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The Hessian of Uj appears here; we denote it by 

in the subsequent. By combining and ()H.4|1 . we find 

{Vsj^q = (det G)Aj ^(VZ\7^)'^G'Vm/ Sf 3^ + 

f 

(det G)Z\t1 J^iVufyCHjGVu^ Sf s„ 

which we rewrite as 

{Vsj^q = (detG)x 

J2 [Aj\VukVGH,GVuk] sl+ 

(3.5) 

J2 [AjiVAj^ydVuk + 2Aj\Vu,yGH,GVuk\ Sk s,+ 

J2 [Aj\\/ukVGHjGVui\ Sk 

to single out the coefficients of different products of s^- 

In the next four propositions, we show that the coefficients of mixed 
products vanish, and that the other coefficients are of the form 

-lA((detG)4.). 

For the proofs of these propositions, we need an elegant representation 
formula for the gradients of the generalized metric coefficients. This formula 
also yields useful identities when extended to the "off-diagonal" coefficients. 

Lemma 3.2 (Representation formula). The gradient of {S/UjY GVuk 
can he represented as 

V{5jkAj) = HjGVuk + HkGVuj + ((Vw^yiV) Vm^ + {{Vuk)^ N)Vuj. 

Thus 

WAj = 2HjGVuj + 2{{Vuj)'^ N)Vuj 

and 

HjGVuk + HkGVuj + {{VujyN)Vuk + ((Vufc)"^iV) Vu^- = 0, k. 
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Proof. By Leibniz' rule, 

^ dqa dqb ^ dqa dq, 



aJ 



•b 



The relation dGah/dq^ = 6aiNb + SbiNa now gives the desired representation 
formula, from which the other two follow immediately. □ 



Proposition 3.3. For all j = 1, . . . , n, 

du 



Aj{VA-^YGVuj + A-^\VujyGHjGVuj = -\Aj^\og{{detG)Aj). 

"''i 



Proof. By writing the left hand side as 

-\{yujyG{2V\ogAj - 2Aj^HjGVuj), 
we can use Lemma (3.21 and p.2|) to further rewrite it as 

- i(VMj)"^G(VlogZ\j- + 2Aj\Vu,y{^GV\ogdetG)Vuj) = 

- liVujyGV {log Aj + logdetG), 

which, by Lemma (3.11 equals the right hand side in the proposition. □ 
Proposition 3.4. For all distinct j, k = 1, . . . ,n, 

Aj\VukVGH,GVuk = -lAk-^\og{{detG)Ak). 

Proof. Lemma and the generalized orthogonality ()1.5p shows that the left 
hand side equals 

- Aj\VukyG[HkGVu, + {{Vu,yN)Vuk + {{Vuky N)Vu,] = 

- Aj' [{Vuk)'^GHkGVu, + Aki^Uj^N] . 

By using Proposition 13. 5| we can further rewrite it as 

- lAj'Ak [(V log Ak^GVu^ + 2{VujyN] = 

- ^A-'AkiVu,)'^ [GV \ogAk + 2N] , 

which, by (IL2j) and Lemma l3?H equals the right hand side in the proposition. 

□ 
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Proposition 3.5. For all distinct j,k = 1, . . . ,n, 

AjiVAjyCVuk + 2A-\VujyGHjGVuk = 0. 

Proof. The left hand side can be rewritten by applying Lemma 13.21 to the 
factor Aj(VA~^)'^ = —Aj^(VAj)'^, after which the result will vanish because 
of the generalized orthogonality p.5|) . □ 

Proposition 3.6. For all distinct j,k,i = 1, . . . ,n, 

{VukVGHjGVue = 0. 

Proof. Lemma f3.2l states that 

HjGWuk + HkGVuj + ((VMj)^iV) Vwfc + {{VukY N)Vuj = 0. 

If this is multiplied from the left by (VuiyG, we have, due to the generalized 
orthogonality ()1.5|) . 

{VueyGHjGVuk + {VuiYGHkGVuj = 0. 

By repeating this process for all other combinations of indices, we get two 
more equations of this form. Together with three trivial symmetry relations 

{VugyGHjGVuk - {VukydHjGVue = 

that arise by transposition of the involved scalars, we get in total six homoge- 
neous algebraic equations for six unknowns of type {VukYGHjGVui. This 
system is uniquely solvable, so the only solution is the zero solution. □ 

The preceding four propositions together with ()3.5|) finally completes the 
proof of (jS3). 

3.3 Stackel form 

In the following two propositions we show that the Hamiltonian ()2.8|1 has the 
Stackel form 

n ^ 

where fk and gk are some functions of one variable only. 
Proposition 3.7. 

(detG)zi. = ^ 
for some functions fk of one variable. 
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The proof follow the proof of Proposition 32 in |H] . 

Proof. We will prove that V[(det G)U'{uj)Aj] = XjVuj for some function Xj. 
This identity implies the proposition, for if it is multiplied by {'Vuk)^G from 
the left, it follows from Lemma ITTl together with the generalized orthogonal- 
ity ()1.5p that {det G)U' (uj) A j only depends on Uj. 

In order to establish the above identity, let us calculate the gradient of 
the identity det(G — ujG) = 0; we get 

d 

(3.6) — det(G - fiG)Vuj + V det(G - /iG) = 0, when /i = Uj. 
For the first term we have immediately 

^ det(G - fiG) = (-l)"(det G)U'{uj), when /x = uj. 

For the second term we use ()1.2|) applied with the pencil — to get 

(3.7) Vdet(G-^G) = 2coi{G - ^G){N ~ ^iN). 
Next, we multiply (jSH) by 2{N^ - X) from the left to find 

(-l)"(det G)U\uj)A.j + 4(A^ - w^iV)"^ cof (G - UjG){N - ujN) = 

(we postpone the technical details to the end of the present proof). We have 
now control over the gradient of {det G)U'{uj) A j; it is proportional to 

V[(iV - UjNy cof(G - UjG){N - ujN)] = 

-^[{N - /iiV)T cof(G - fiG){N - iiN)]Vuj+ 

V[(iV - jjNy cof(G - fiG){N - 12N)], when /j = uj. 
In fact, the last term here is also proportional to Vuj, since 
V[(iV - fiNy cof (G - fiG){N - fiN)] = 2{a - /la) cof (G - iiG){N - fiN) 
and 

d 



2 cof (G - iJ,G){N - ij,N) = — det(G - iJ.G)Vuj, when n = u 
as follows from ()3.7j] and ()3.(i|l 
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To complete the proof, we need to show two identities. The first one is 
2(A^"^ - N'^X)Vuj = Ay From apphed with both and N and from 
fTTj) we find 

2{N^ - N^X)Vuj = [(VlogdetG)"^G - (V logdet G)"^G'X]Vmj- = 

(VlogdetX)"^GVM^- = {VtiXyx-^GVuj = (VwjO'^G'VtrX, 

which by Lemma f3. II equals Aj. 

The second identity we need to show is 

{N^ - N^X) cof{G - UjG){N - u^N) = 

{N - UjXy cof(G - UjG){N - UjN), 

which is immediate once we have observed that (X — fil) cof{G — fiG) = 
det(X — fil) cof G vanishes when fi = uj. Indeed, 

{N^ - N^X) cof(G - i2G){N - fiN) = 
{N - fiXy cof (G - nG) {N - fiN) -N^{X - fil) cof (G - fiG) {N - fiN) , 

which in turn equals {N — fiN)'^ cof(G — fiG){N — fiN) when fi = Uj. □ 
Proposition 3.8. 

(3.8) ^ = g_J_,,(„.), 

for some functions gi of one variable. 

Proof. The idea is to solve the Frobenius compatibility conditions for Equa- 
tion ()1.4j) satisfied by any pair of quasi-potentials belonging to the same 
cofactor pair system. In terms of u coordinates, this equation can be written 
as XJduW = {det X)JduW, or 

UduW = {detl{)duW 

with U = J^^XJ as in p.6|) . Since lA is diagonal, we find 

= ( TT"« ) k = l,...,n. 

For W to exist (we assume that it is in C°°), we require its mixed second 
order partial derivatives with respect to k and i to be equal. This gives 

dW d'^W dW d'^W , , 

^ + ^^ r. ^ =^ + ^fc a a ' fc,£ = l,...,n. 

ouk oukOUi oue ouiOUk 

These equations can be written in the form d'^{uk — ue)W/dukdue = 0, which 
according to CoroUarv IB. 21 have the indicated solution. □ 
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3.4 Quadratic first integrals 

We shall now show that the functions 



n 



i=l 



U' U. 



with fi from Proposition 13.71 and Qi from Proposition are first inte- 
grals for the Hamiltonian system governed by the Hamiltonian ()2.8p . that 
is, H = E^"~^\ This fact can be expressed in terms of the Poisson bracket 
as {H,E^''^ = 0. Since H has the same form as all other functions E^''\ it 
requires no extra effort to show that {E^^\ E^} = for all k, i. This means 
precisely that all first integrals E^''^ are in involution. 

Proposition 3.9. The functions E^''^ are first integrals in involution. 

Proof. A straightforward calculation shows that 



i=l 



a=l 



U' 



fi{Ui) (Ti{Ua) - T:— ———fi[Ui) ak[Ua) 



dUa \U'{Ui) J dUa \U'{Ui) 

-gi[ui) ai[ua) - ^ — TT7r~^9i[ui) ak{Ua) 



+ 



dUa \U'{Ui) J dUa \U'{Ui 

This is identically zero since the expressions within square brackets vanish. 
The only property of the functions fi and gi needed to see this is that they 
only depend on one variable each. 

If a = z, we have dak{ui) / dua = dai{ui)/dua = 0, so that it is possible to 
take these symmetric polynomials outside the respective derivative, showing 
that the two terms cancel. 

If a 7^ 2, we can factor out fi{ui) or gi{ui) and thus have to show that 

d [ (Tk{ui)\ . d f ae{ui)' 



dUa\U'{Ui)j dUa\U'{Ui 

Since d[l/U'{ui)) /dua = 1/ [U'{ui){ui — it is equivalent to show that 
d 

{Ui - Ua)^(yk{Ui) Cri{Ua) + Ok{Ui) Cr^(Ua)- 

d 

OUa 
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From ()A.3|) and ()A.4|) follows that the left hand side can be written as 



{ui - Ua)o-k-i{uiai>^i + ai) + {uaOk-i + ak){ui(Ji-i + a^)- 

{Ui - U„)o-£_i(Ui(Tfe_i + O-fe) - {UaCft^l + 0-^)(Ui(Tfe_i + (Tfc) 

in terms of elementary symmetric polynomials depending on all mi, . . . , u 
except Ui and Ua- This expression is easily seen to vanish. □ 



n 



The first integrals E'^^'^ are functionally independent as functions on i?^" = 
{(m,s)}. The standard test for this is that the differentials dE^^^ should be 
linearly independent. Here it is simple to verify this fact, since the matrix 
{dE^^'^ / dsi) can be written as a product of two non-singular matrices: the 
inverse Vandermonde matrix ()A.5|) and diag(/i, . . . , /„). 

3.5 Solution by quadratures 

We shall now use our previous results to obtain the solution of the cofactor 
pair system ()2.4j) by quadratures. 

We begin by noting that the first integrals separate in the following sense. 
Denote by the constant value of E^^\ and consider the equations 

(3.9) E^^\u,s) = ak, fc = 0,...,n-l. 

These equations can be put into matrix form as 

/ (To(mi) o-o(m„) \ 



U'{ui) U'{Un) 

Cr„_i(Mi) an-l{Un) 



\fi{ui) sl + gi{ui) 

J/n(Mn) sl + QniUn) 



H 


'Oin-l\ 











V f/'(ni) f/'K) / 

from which it is clear that they are equivalent to 

= (-1)"+^ 

\fn{Un) + gn{Un)) ^ ... {-Unf ) \ 

Hence ()3.9|) is equivalent to the system 

n-l 

(3.10) \fk{uk)sl + gk{uk) = {-lT^^Y.^-UkYai, A; = l,...,n, 
consisting of relations involving only the A;:th coordinate pair (m^, Sk)- 
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We then define a function S{u) = ^ Sk{uk) by requiring the functions Sk 
to be solutions to their respective "separation equations" 



2 



n-1 



i=0 



The function 5* so constructed will thus be a simultaneous separated solution 
to the Hamilton- Jacobi equations 



In particular this holds for the Hamiltonian H = E^^~^\ that is, ()2.8|) . 
By means of the classical Hamilton-Jacobi method we can therefore obtain 
the solution (M(r),s(r)) of the Hamiltonian system corresponding to H by 
quadratures. By transforming this solution to variables, we also have 

the solution (g(r),p(r)) of the system ()2.5p . 

Finally, to get the solution of the cofactor pair system ()2.4|) . we need to 
find the function r(t) from the differential equation ()2.9jl . which requires one 
further quadrature. The solution for the cofactor pair system can then be 
obtained as (g(r(t)),p(r(t))). 

3.6 Conclusion 

In Section 13.11 and Section 13.21 we established a canonical Hamiltonian form 
for the system ()2.5|1 related to the cofactor pair system ()2.4|1 . In Section 1331 
we showed that the corresponding Hamiltonian has Stackel form, which made 
it possible to construct n quadratic first integrals in Section 13.41 In Sec- 
tion 13.51 we used the Stackel form to separate variables in the corresponding 
Hamilton-Jacobi equation, providing the solution to the related system. The 
solution to the cofactor pair system was then obtained by rescaling the time. 

This completes the proof of Theorem 12.11 

4 Separation coordinates for generic cofactor 
pair systems 

In the separability theorem we defined the separation coordinates Uk for 
the cofactor pair system cPq/dt^ = — (cof G')"^VH^ = — (cof G')"^VW^ as 
roots of the polynomial det(G' — fxG), provided that the functions Uk{q) are 
functionally independent. Equivalently it can be said that Uk are eigenvalues 



E^''\u,duS) =ak, k = 0,...,n 



- 1. 



of X = G'^G, since det G 0. 
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In order to better explain the nature of these separation coordinates, we 
shall derive explicit formulas for the coordinate surfaces, which are "non- 
confocal quadrics". These formulas show how these new separation coordi- 
nates relate to the well known elliptic and parabolic coordinates that fre- 
quently occur in the classical separability theory. 

4.1 Classical elliptic and parabolic coordinates 

When G = I, the cofactor pair system can be thought of as a classical 
potential Newton system (Pq/dt^ = —VW admitting an extra quadratic first 
integral E = ^q'^ {cof G)q + W in addition to the energy integral. In the 
generic case, the separation coordinates Uk are then eigenvalues of G, which 
define elliptic or parabolic coordinates as explained by the following three 
propositions [7j. 

The first proposition shows that the matrix G can be given a particularly 
simple standard form by a change of Euclidean reference frame q ^ Sq + v 
(such an affine transformation is called Euclidean if is an x n orthogonal 
matrix, S'^S = I, with det S = 1, and v G i?" is a translation vector). Under 
this transformation the quadratic forms 

q^ coi[G{q)]q and q'^q 

go into 

ci^ coils'^ G{Sq + v)S]q and q'^q 

so that we still have a potential Newton system with an extra quadratic 
first integral. One can then give S'^G{Sq + v)S one of the stated standard 
forms by a suitable choice of 5* and v, as in the proofs of Theorem 14.51 and 
Theorem 14.101 below. 

Proposition 4.1 (Standard form). Let G{q) = aqq^ + I3q^ + q(3^ + , with 
a ^ R and (3 G -R" not both zero, he an elliptic coordinates matrix. There is 
a Euclidean change of reference frame q ^ Sq + v , which gives G one of the 
following standard forms: if a ^ then 

G{q) = -qq^ + diag(Ai, . . . , A„), 

and if a = but (3^0 then 

G{q) = e„g"^ + qe^ + diag(Ai, . . . , A„_i, 0), = (0, . . . , 0, 1)"^. 

The following two propositions can easily be proved by using a Weinstein- 
Aronszajn formula, as in the proofs of Theorem 14 . 51 and Theorem 14 . 1 ( )l below. 
These propositions justify the name elliptic coordinates matrix for G. 
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Proposition 4.2 (Elliptic coordinates). // 

G'(g) = -qq^ + diag(Ai, . . . , A„), 
then the eigenvalues Ui{q), . . . , Un{q) of G satisfy the rational equation 

^2 



A. 

i=i ' fc=i 



(4.1) 1 + = n(/^ - ^^■) / n('" - ^'^)' 
i=i ^ * 



which, if Xi < ■ ■ ■ < Xn, is the defining equation for elliptic coordinates uj 
with parameters Xk- 

Proposition 4.3 (Parabolic coordinates). // 

G{q) = (^nq^ + qeZ + diag(Ai, . . . , A„_i, 0), e„ = (0, . . . , 0, 1)"^, 

then the eigenvalues ui{q), . . . , Un{q) of G satisfy the rational equation 

n— 1 2 " / '^^1 

(4.2) 2g„ - /i + 5^ = - n('" - 

i=i ^ * j=i ' k=i 

which, if Xi < ■ ■ ■ < Xn-i, is the defining equation for parabolic coordinates Uj 
with parameters Xk- 

Equation ()4.ip define coordinate functions Uk{q) as zeros of the rational 
function (of /i) in the left hand side. It is easy to see that this function has 
n distinct real roots ui{q), . . . , M„(g) such that 

■Ui < Ai < ^2 < A2 < ■ ■ • < Mn < A„ 

for each n-tuple (gi, . . . , g„) of non-zero reals qk- All these roots Uk clearly 
satisfy 



E 



9? 



A,- — u 



i=l 

from which it follows that the hypersurfaces u{q) = c of constant value of 
u{q) are confocal quadrics (quadratic surfaces) in R"'. Their character is 
decided by the value of c: for c < Ai, the equation determines a family of 
ellipsoids, for Ai < c < A2, a family of one-sheeted hyperboloids, and so on, 
over various families of hyperboloids, until A„_i < c < A„ (when c > A„, 
there are no real solutions for q). The functions Uk{q) are then called elliptic 
coordinates. 
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Similarly, Equation ()4.2|) define parabolic coordinates Uk{q) satisfying 



Ui < Ai < ^2 < A2 < ■ ■ ■ < A, 



'n-l < U. 



'71 



and 



71-1 




u + 



2g, 



i=l 



which also define a family of confocal quadrics. 

In both cases, the coordinate surfaces are orthogonal to each other as 
follows from the generalized orthogonality condition {'Vuj)'^G'Vuk = with 
G = I. That is, the elliptic and parabolic coordinate systems are orthogonal 
coordinate systems. 

In general, for arbitrary cofactor pair systems G ^ I, and the separation 
coordinates are non-orthogonal with respect to the standard scalar product. 
However, they are always orthogonal with respect to G{q). 

4.2 Cofactor-elliptic coordinates 

We shall now derive rational equations similar to 1)4.11) and ()4.2|) for coordi- 
nates defined by two arbitrary elliptic coordinates matrices G and G. To this 
end, we shall use a formula for the Weinstein-Aronszajn determinant 

det{I + AR), 

where A and R are n x n matrices. If R is non-singular and A is of low 
rank m, this determinant is effectively of dimension mxm, which can be seen 
by changing to a suitable basis. However, the factorised representation AR 
is not important here, so we will instead consider det(J + Ai), where Ai is 
of low rank. (The general case then follows by observing that Ai = AR has 
the same rank as A.) 

Proposition 4.4 (Weinstein-Aronszajn). Let A4 be an n x n matrix of 
rank m < n. Given an orthonormal set {/i, . . . , fm} of vectors in R" that 
spans the range of M., the Weinstein-Aronszajn determinant can he computed 
as an m X m determinant: 



det{I + M) = det{S,, + fjMfk) 
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In particular, for m = 1 



det(/ + A^) = l + fjMfu 
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and for m 



det(/ + M) 



1 + fjMh fjMh 
fjMh 1 + fjMf2 



Proof. Extend the given set {/i, . . . , /„} to an orthonormal basis {/i, . . . , fm, 
fm+i, ■ ■ ■ , fn} for R"'. We represent ^A hj a. sum of m outer products as 
= Y^a=i fa fl with fa = M'^ fa- In the standard basis {ei, . . . , e„}, the 
corresponding matrix is then F'^MF = Yl^=i ej, where F = (/i, . . . , fn) 
and Ca = F'^ fa- Thus in F'^ AiF all elements in rows m + 1, . . . , n vanish. 
This implies the block form 



det(/ + M) = det(/ + F^MF) 



h + fJMfj 







5, 



which shows that expansion of this determinant according to the last row 
results in a new, lower-dimensional determinant of the same form. This 
process can thus be repeated to eventually find an m x m determinant. □ 

We now apply this proposition to the case when at least one of G and G 
is quadratic in g, that is, at least one of a or a is non-zero. 

Theorem 4.5 (Cofactor-elliptic coordinates). Let G{q) = aqq^ + qP'^ + 
(3q^ + 7 and G{q) = aqq^ + q[3^ + [3q^ + 7, with not both a and a zero, he 
two elliptic coordinates matrices. Define a vector B = B{^) in and a 
symmetric n x n matrix F = F{^) as 



a 



r(/i) 



7-/^7 



a 



Choose an orthogonal matrix S = S{^) so that S'^{r — BB'^)S is diagonal 
with eigenvalues Ai = A„_i = A„_i(;u). Define a fi-dependent 

Euclidean change of reference frame by Q{q) = S'^{q + B). 

If all Aj(/i) are non-zero, the roots Ui{q), . . . , Un{q) of the polynomial 
det(G' — fiG) satisfy the rational equation 



Qijq;!^)^ _ det G{q) 
Aj(/i) (a — /ia)" 



-/i) / JjAfc(/i). 



(4.3) 

— A. I /A I \ (y — ii,(y f 

j=l ' k=l 

Proof. In order to use Proposition 14.41 we rewrite the pencil G — fiG as 

G - fiG = {a - ^ia)qq^ + {(3 - ^if3)q^ + q{(3 - ^I3Y + 7 - /i7 = 
[a - ^la) [{q + B{^^)) [q + i?(/x))"' + r(/x) - i?(/i) B{^Y] = 

(a - ^ia) [A{q; /i) + T(/i)] = (a - /ia) [/ + A{q; /i) T(/i)-i] T(/i), 
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where A = {q + B){q + and T = T - BB'^. Using that G - /iG = 
G(X -/iJ), we find 



(4.4) 



det(/ + AT'^) 



detG detjX -fil) 
(a — /ia)" det T 



The left hand side can easily be calculated by invoking Prop osit ion 14 . 41 f or the 
rank 1 matrix AT~'^, and by diagonalizing the symmetric matrix T. Indeed, 
by setting /i = {q + B)/\q + B\, and by choosing an orthogonal matrix S 
such that S'^TS = A = diag(Ai, . . . , A„), we have 



det(/ + AT-^) = 1 + fjAT'^fi = 

l + {q + BYT-\q + B) = I + {S'^ {q + B))" A'^S'^ {q + B)) . 



As is indicated in the following remark, the coordinates that satisfy ()4.3p 
generalize the classical elliptic coordinates. Since these new coordinates nat- 
urally arise in the theory of cofactor pair systems, we call them "cofactor- 
elliptic" coordinates. 

Remark 4.6 (Elliptic coordinates). If Theorem 14.51 is applied with 



then i? = and F = diag(/i — Ai, . . . , fi — Xn), and we recover Proposition lOl 

Example 4.7 (Cofactor-elliptic coordinates in R^). This example has 
previously been studied by H. Lundmark, who calculated the relevant ex- 
pression directly from the pencil G — fiG by completing squares. 
Suppose that 



where Ai, A2 and e are parameters. The case of one quadratic and one linear 
elliptic coordinates matrix is in fact the most general, since if both G and 
G were quadratic, one could consider the pair G and G — {a/a)G instead. 
When £ = 0, we have the situation in Remark I4.(i[ 

Following the notation of Theorem 14. 5| we have B = (e/x, 0)""^ and F = 
diag((l — e)fi — Ai, /i — A2) , so that 



Equation ()4.H|1 now follows from ()4.4j) . 



□ 



G{q) 



qq + diag(Ai, . . . , Art) and G = I, 




and 




F — BB^ = diag((l — £)fj, — Ai — e^/i^, fi — X2) and S = F 
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From the theorem we get the rational equation 



(1 — e)^ — Ai — e^/i^ /i — A2 



((1 - £:);U - Ai - eV) - -^2) ' 

describing the relation between the separation coordinates (mi,M2) and the 
Cartesian coordinates (gi,g2)- The coordinate curves of constant = u are 
quadrics given by 

(A ^) (gi + ql ^ 

^ ' \i-{l-e)u + e^u^ \2-u 

This relation is more complicated to comprehend than the corresponding one 
for elliptic coordinates. It is the signs of the denominators that determine 
the kind of quadrics involved. Here we consider two cases: e = and e = 1. 
When £ = 0, ()4.5p specializes to 

(4.6) + = 1. 

\l — U A2 — U 

We assume that Ai < A2 (there is no qualitative difference in assuming 
Ai > A2) and study the signs of the denominators as u varies. For m < Ai, 
we have always (+, +) defining a family of confocal ellipses; for Ai < m < A2, 
we have (— , +) defining a family of confocal hyperbolas; finally for u > A2, 
we have (— , — ) for which there are no real solutions (gi,g2) to ()4.(i|l . These 
two families of curves define 1-1 mappings of the region {(ui, U2) G i?^; ui < 
Ai < U2 < A2} onto each quadrant of the gig2-plane; therefore we have a 
coordinate system in the whole plane except on the qi- and g2-axes. 
When e = 1, ()4.5p specializes to 

(4.7) 



Ai + A2 — M 

and the picture is more complicated. First of all, if Ai > 0, then Ai + 
> 0, and we have (+, +) if m < A2 and (+, — ) if m > A2. Also in this 
case the equation define a family of ellipses and a family of hyperbolas, but 
they are not confocal since the curves are being translated along the gi-axis. 
The curves do not fill the entire gig2-plane either, and we can only define a 
coordinate system in certain regions of the plane. Figure ^ illustrates this 
situation. 
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-10 -5 5 10 

Figure 1: Solutions to ()4.7|) with Ai = 5, A2 = and — 5 < m < 5. 



Case 



Signs and curve types along the u axis 



A2 G (-00, -A^; 



(+,+) 



Ell 



A, 



(+,-) 



Hyp 



-A^ 

Void 



(+,-) 



Hyp 



\2e{-\l,X\) 



> —\l > A2 > Ai 



Ell 



Hyp 



Void 



Hyp 



A2 G (A];, +00) 



( + : + ) (-. + ) ( + : + ) ^ (+.-) , 
^ —\l ^ \l > A2 > 



Ell 



Hyp 



Ell 



Hyp 



Table 1: Different solution types to ()4.7p with Ai < 0. Notation A^ = 
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-10 -5 5 10 

Figure 2: Solutions to ()4.7|1 with Ai = —5, A2 = —3 and — 5 < u < 5. 
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-10 -5 5 10 

Figure 4: Solutions to ()4.7|) with Ai = —5, A2 = 3 and —5 < n < 5. 



The picture is similar when Ai < 0, except that three or four families 
of curves are defined as is indicated by the combination of signs in Ta- 
ble E Three different cases occur depending on the relation between A2 
and ±\/—\i. They are illustrated in Figures El El and IH 

Example 4.8 (Cofactor-elliptic coordinates in R^). We now extend 
Example 14.71 in the case e = 1 to three dimensions. We take 

/-qf + Xi -qiq2 -QiQ-i \ _ /2gi q2 

G{q) = -qm -qh + A2 -q2qz and ^(g) = \ q2 1 

V -?i?3 -^2^2 -gi + Aa/ \qz 

where Ai, A2, A3 are parameters. As before, we find B = (/i, 0,0)""" and 
r = diag(— Ai, — A2 + /i, —A3) so that 

r — BB^ = diag(— Ai — /x^, — A2 + /x, —A3) and S = I. 

We get the rational equation 

1 , , 4 4 _^2 {Ui- fl){u2- fi){u3- 

+ \ ,,2 + \ , , + \ ~ VS" 



-Ai-/i2 -A2 + /X -A3 "'(-Ai-/i2)(-A2 + /x)(-A3)' 

describing the relation between the separation coordinates {ui,U2,U3) and 
the Cartesian coordinates {qi, q2, q'3). The curves of constant n = u satisfy 

(gi + uf ^ ql ^ ^ _ ^ 



Ai + A2 — M A 



3 
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If A3 > 0, there are always real solutions (q'l, ^'2, Q'3) to this equation for any 
given u. The character of the curves are determined by the relation between 
Ai, A2 and u, as explained in Example 14.71 

4.3 Cofactor-parabolic coordinates 

We shall now investigate the separation coordinates in the case when none 
of G and G are quadratic, but at least one is linear, in q. That is, we assume 
here that a = a = and that one of f3 and [3 is non-zero. As in the case of 
Theorem 14. 5[ we start with Proposition 14.41 In order to formulate and prove 
a similar theorem, we need the following fact about "partial diagonalisation" 
of symmetric matrices. 

Proposition 4.9 (Partial diagonalisation). Let F be a symmetric n x n 
matrix, and let {Bm+i, ■ ■ ■ , -Bn} be an orthonormal set of vectors in R^. There 
exists an orthogonal matrix S with columns Sj = Bj for j = m + 1, . . . ,n, 
such that S'^rS attains the block form 



where D is a diagonal m x m matrix, Ci is an m x {n — m) matrix and C2 
is a symmetric {n — m) x {n — m) matrix. In particular, if m = n — 1, there 
is a vector c such that 

S'^rS = diag(Ai, . . . , A„_i, 0) + e„c"^ + eel, e„ = (0, . . . , 0, 1)"^. 

Proof. Take any orthogonal matrix S with columns Sj = Bj for j = m + 
1, . . . , n. Then 



for some matrices D, Ci and C2. Choose an orthogonal m x m matrix P 
such that P'^DP = D and set 






Then S has the stated last columns, and 




where Ci = P'^Ci and C2 = C2- 



□ 
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Theorem 4.10 (Cofactor-parabolic coordinates). Let G{q) = qjS^ + 
/3g"'" + 7 and G{q) = qP'^ + Pq'^ + with not both (3 and (3 zero, he two elliptic 
coordinates matrices. Define a vector B = B{^) in i?" and a symmetric nxn 
matrix F = r{^) as 



1/3 - ^^p\ 1/3 - ^ip\ 

Take S = S{fi), Ai = Ai(/i), . . . , A„_i = A„_i(/i) and c = c{fi) from Propo- 
sition HiHl applied with F and {B}. Define a fi-dependent Euclidean change 
of reference frame by Q{q) = S'^q + c. 

If all Aj(/i) are non-zero, the roots Ui{q), . . . , Un{q) of the polynomial 
det{G — fiG) satisfy the rational equation 

,,8) 2g„(, ,) - £ . . n(., - ,)/ n M,), 

Proof. The proof is similar to that of Theorem 14. 5t but more technicaL We 
rewrite the pencil G — ^G as 

G-fiG={P- fiP)q'' + q{(3 - fipy + 7 - /i7 = 

1/3 - /i/3| + Vy + {q + V)B{^y + r(^) - B{^) -V i3(/x)"^) = 

1/3 - /x^l {A{q- fi) + Tifi)) = 1/3 - /x/31 (/ + Aiq; ^) r(/i)-i)r(/i), 

where A = B{q + + {q + V)B'^ = BqJ, + qyB'^ , which defines qy, and 
where T = F — BV'^ — VB'^ and V is an as yet undetermined vector. We 
get 

/. N -, In det G det (X — uJ) 

4.9 det / + AT-^) = . \ ^ ' . 

^ ' ^ ^ |/3-/i/3|" detr 

The left hand side can be calculated by invoking Proposition 14.41 for the 
rank 2 matrix AT~^ by setting fi = B and /2 = {qv — {QvB)B) /a where 
'^^ = QvQv — iOvBy. A straightforward calculation gives 



1 + fjAT-'f, fjAT'^f2 
fjAT-^f, 1 + fjAT~^f2 



[I + qlT-'Bf - {q^T-'qy){B^T-'B) 



This is the point were it should be apparent that an ordinary diagonal- 
isation of T does not help; the problem is the product T~^B, which would 
attain the form S S'^ B , which is not very useful unless we can simplify the 
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expression S^B. However, Proposition 14.91 gives us an orthogonal matrix 5* 
having B as last column, and therefore S'^B = e„, which proves helpful. 

The matrix S provides the decomposition S'^^FS = diag(Ai, . . . , A„_i, 0) + 
Cnc""" + cej[, which implies 

S'^TS = diag(Ai, . . . , A„_i, 0) + e„(c - S^V^ + (c - S''V)el. 

We now want this to be diagonal with non-zero diagonal elements, since 
we already have assumed that exists. Therefore we choose V so that 
c — S'^V = |e„, that is, we choose V = Sc— \B. Thus 

5TT5' = ^ = diag(Ai,...,A„_i,l), 

and the above expression for det(/ + AT~^) becomes, with Q = S'^qv, 

""^ n2 

(1 + Q^A-'erf - iQ^A-'Q)ielA-'en) = 1 + 2g„ - V ^. 

Changing to the notation in the formulation of the theorem, Qi = Qi for 
i = 1, . . . , n - 1 and 1 + 2Q„ = 2Q„, which shows that gives □ 

As in the elliptic case, the coordinates defined by ()4.8|1 will be called 
"cofactor-parabolic." 

Remark 4.11 (Parabolic coordinates). If Theorem 14. 101 is applied with 

^(9) = en?""" + ge^ + diag(Ai, . . . , A„_i,0) and G = /, 

then B = Cn and F = diag(Ai — /i, . . . , A„_i — /i, — /i), so that S = I and 
c = — i/ie^. This gives Q = q — |/ue„, that is, Qi = qi for i = 1, ... ,n — 1 
and 2Qn = Iq^ — /i. Proposition 14.31 is thus recovered. 

Example 4.12 (Cofactor-parabolic coordinates in I^). Suppose that 

where A is a parameter. Following the notation of Theorem 14.101 we have 
B = (-/i, l)T/A^(/i) and F = diag(A, 0)/iV(/x), where iV(/i) = ^1 + /i^. In 
Proposition 14.91 we can take 
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so that 

cT p c — ( ^ "Z^"^) 

from which Ai = \/N{fiY and c = (— /uA, |/i^A)"'"/A^(/i)^. The theorem gives 
the rational equation 

iV(/x) +iV(/x)3 A V ^(/^) Ar(/i)3j 

N{^iY ' A 

For a fixed fi = u, we thus have 

g g2 - ^gi M^A _ iV(n)^ / gi + wga _ V _ 

Ar(M) ^ Ar(u)3 A V N{u) N{u)^ / ~ 

which after multiphcation by N{u) and expansion of the square simphfies to 



(4.10) 2g, - i^l±!^ = 0. 

A 

This equation defines a fan of parabolas centered at the origin, filling the 
entire upper half plane when A > 0, and similarly in the lower half plane 
when A < 0. This is illustrated in Figure El 

Equation ()4.1()|1 allows for a global definition of the coordinates, as is 
easily seen by solving the equation for u. Indeed, for a fixed A we can take 

-qi + V2Ag2 , , s -qi- 

uAq) = and U2{q) = 

g2 g2 

everywhere in the region {{qi, ^2) ^ R^', ^12 > 0}. 



4.4 The functional independence assumption 

We complete this study of the separation coordinates Uk by giving a formula 
that can be used to investigate whether the functions Uk{q) are functionally 
independent without having to calculate them explicitly. Recall that they 
are functionally independent if det J 7^ 0, where J = (Vmi, . . . , V-u„) is 
the previously defined Jacobian. The formula is formulated in terms of the 
coefficients in the polynomial det(X — ///), but if only G and G are given 
it is not necessary to calculate X = G~^G, since det(X — fil) = det(G' — 
/iG)/detG. 
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Figure 5: Solutions to ()4.10|) with A = 1 and — 5 < m < 5. 



Proposition 4.13 (A formula for the Jacobian J). Let a^, . . . , an he 

the coefficients in the characteristic polynomial of X defined by 

n 

det(X - nI) = Y^ an-k 

fc=0 

and let Uk he the roots of this polynomial. Then 

(4.11) det(Vai, . . . , Van) = det(Vui, . . . , Vu„) ■ - Uk). 

j<k 

Proof. Since = ak{u) for all k, the chain rule and ()A.3jl imply 

sr^ dak(u) \ ^ 

j 3 

These relations can be put into matrix form as 

/(Jo(mi) ... cr„„i(Mi)\ 
(Vai, . . . , Va„) = (Vmi, . . . , Vu„) ■ j . 

\a-o(M„) ... (Tn-liUn)/ 

This implies (ll.llj) . for the matrix [(Jjiuk)) is the transpose of the matrix 
()A.7j) . which according to Appendix 1X1 has determinant nj<fc('"i ~ ^fc)- ^ 

Remark 4.14 (An a priori functional independence check). It is pos- 
sible to use the formula also to check if the quadratic first integrals defined by 
a cofactor pair system are functionally independent. Indeed, in view of the 
separability theorem, they are functionally independent if the functions Uk{q) 
are, and the formula shows that it suffices to check the coefficients ak{q). 
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A Symmetric polynomials, etc. 

In this paper we frequently encounter two types of elementary symmetric 
polynomials. These are cri{u), which are defined by 



(A.l) l[{z + Ui) = J2^'"~'^^ 

i=l 

and (Ti{uj), which are defined by 



i=l i=0 



n.-l 



The symbol Uj indicates that <Ji{uj) depends on all mi, . . . , except Uj. We 
also extend these definitions to include a_i{u) = and a_i{uj) = 0. 

There are two sets of identities that provide a connection between these 
polynomials. The first is 

d 

(A. 3) —ai{u)=cri_i{uj), z = 0,...,n, 

which follows by comparing the derivative of ()A.1|) with ()A.2|) . The second 
is 

(A.4) ai{u) = Uj(Ti^i{uj) + ai{uj), i = 0,...,n, 

which follows by expanding the product in the left hand side of ()A.1|) . except 
for the factor {z + Uj), using ()A.2|) . 

Closely connected with these polynomials are the products 



dz 



W{Z - Ui) = Y\_iuk - Ui). 

Z=Uk i=i i^k 



They arise in the inverse 

(yQ{ui)/U'{ui) ... ao(M„)/f/'(M„) 

^cr„_l(Ml)/f/'(Ml) ... CT„_l(M„)/f/'(M„) 

of the Vandermonde matrix 



(A.5) V-^ 



(A.6) V=(-l 



-uiY-^ ... (-ui)O' 
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By using ()A.2|) . it is easy to check that the two matrices are each others 
inverses: the element in position (j, k) in the matrix product V ■ is 



i=l 

f-r 



n+l /I \2n 



jk- 



We recall the well known fact det V = nj<A;('"j " ""fc)- We also need the 
determinant of the matrix 

(O-o(^il) • • • (^o{Un) \ 
= V-i-diag(f/'(Mi),...,[/'K)), 
an-iiui) ... ct„_i(m„)/ 

which consequently can be calculated as 

n , 

j=l j<k j<k 

B The Levi-Civita equations 

In the proof of Proposition l3.8| we need an explicit expression for the solution 
to the system d'^{ui — Uj)V /duiduj = of PDEs for V . The solution is 
obtained by considering a special case of the following lemma, which deals 
with a system of PDEs first derived by T. Levi-Civita. 

In the formulation of the lemma, we use the concept of a Stackel matrix. A 
Stackel matrix is a non-singular matrix ^ = det^ 7^ 0, whose i:th 

row only depends on Wj, like, for instance, the Vandermonde matrix ()A.6j) . 

Lemma B.l. Suppose that $ is a Stackel matrix with inverse \1/ = [ipij) 
having non-vanishing elements in the first row. The solution of the system 
of (l) PDEs 

d'^V d log dV d log i^udV ii ^ ]) 

duiduj dui duj duj dui 

is 

(B.2) V{u) = J2^ikfk{uk) 

where fi are some functions of one variable. 
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Proof. We will use the Levi-Civita separability condition j^l C] to find the 
solution of ()B.1|) . This condition states that a Hamiltonian H{u,s) on R^^ 
with coordinates (m, s) is separable, that is, its corresponding Hamilton- 
Jacobi equation H{u, duS) = E admits a complete separated solution S = 
Yl Sk{uk), if and only if the equations 

dijH d'H &H - diH diH H + d'^ H diH djH - diH d'H djH = 

are satisfied identically with respect to (m, s) for all distinct i, j = 1, . . . , n. 
Here we write di for d/dui and for d/dsi. When this condition is applied 
to the Hamiltonian H = ^2 '^ikSk — V, we find 

^(^li ^ij dijipik - i^ii ditjjij djipik - 'ipij djipii di^ik)sk- 

{ipii ipij dijV - ipii diipij djV - ipij djipii diV) = 0. 



The vanishing of the terms involving V is equivalent to (jB.ljl . Similarly, the 
coefficient of Sk vanishes if and only if ipik solves (jB.ljl : below we will show 
that this is in fact so. Thus, if V satisfies fIB.lj) . there is a separated solution 
S = Sk{uk) to the equation ^-^ifc dkS — V = 0, which implies (jB.2jl with 
fk{uk) = dkSkiuk). 

Claim: ipik solves (jB.ljl . To prove this, we introduce the notation ip = 



det$ and ipM ~ i^iki^ji ~ i^ui^jk- We note that the product (pipij is indepen- 
dent of Uj since it equals (—I)*"'"-' times the determinant obtained by deleting 
row j and column i in ip. Likewise, the product v^V^h independent of Uk 
and U£ since it equals (— 1)*+^'+^+^ times the determinant obtained by deleting 
rows k, i and columns i, j in (p In the same spirit, we state the identities 



1 



/ p^ _ jjik ^ _J_ f _ p^ 



which imply that 



dij — - — - — = and dij — 



respectively. On the other hand, by Leibniz' rule. 



dijipik + i^u dj [ — j diipik + ipij di y—j dji^ik 



The expression within square brackets thus has to vanish, which proves the 
claim. □ 
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It is easy to check that ()B.2|) satisfies ()B.1|) with arbitrary functions 
fk{uk)- One would therefore hke to draw the conclusion that the general 
solution of (jB.ll) is ()B.2j) with arbitrary functions fk{uk), but it seems very 
hard to show this. From the above proof also follows that (jB.ll) can be put 
into the suggestive form 

duiduj \(p^pu^pij J 

but this does not seem to help in finding the general solution to the equations 
by direct methods. 

However, by applying the lemma with ^ being the Vandermonde matrix 
()A.6jl . we establish the following. 

Corollary B.2. The system d'^{ui — Uj)V / duiduj = has the solution V = 
fk{'^k)/U'{uk) where fk are some functions of one variable. 
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